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ABSTRACT 

Lithium abundances in a selected sample of halo stars have been revised by using 
the new accurate IRFM effective temperatures by Alonso, Arribas & Martinez-Roger 
(1996a). From 41 plateau stars (T e s > 5700 and [Fe/H] < -1.5) we found no evidence for 
intrinsic dispersion, a tiny trend with T e ff and no trend with [Fe/H]. The trend with 
the T e g is fully consistent with the standard Li isochrones of Deliyannis, Demarque 
& Kawaler (1990) implying a primordial value for Li of A(Li)= 2.238 ± 0.012 1(T ± 
0.05 S j, s . The present results argue against any kind of depletion predicted by diffusion, 
rotational mixing or stellar winds. Therefore the Li observed in Pop II stars provides 
a direct and reliable estimate of the baryonic density that can rival other baryonic 
indicators such as the deuterium in high redshift systems. The present upwards revision 
of primordial Li in the framework of SBBN gives at la two solutions for the baryonic 
density: n B h 2 = 0.0062±g;ggi? or Vl B h 2 = 0.0146±°;°°:g. 

Key words: Stars: abundances - Stars: Population II - Stars: fundamental parame- 
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1 INTRODUCTION 

In a seminal paper Spite & Spite (1982) found that the warm Pop II dwarfs (T e a > 5700 K) share a unique Li abundance, at 
variance with their Pop I counterparts, where a large spread in Li abundances is observed. Since no sign of stellar depletion 
was evident on the first dozen stars the Spite plateau was readily interpreted as a signature of pristine Li abundance directly 
related to primordial nucleosynthesis. The original finding has been later confirmed by a number of investigations (Spite, 
Maillard & Spite 1984, Spite & Spite 1986, Hobbs & Duncan 1987, Rebolo, Molaro & Beckmann 1988). The tight connection 
between Pop II Li abundance and BBN is of far reaching cosmological importance since in the standard BBN the primordial 
Li abundance is a function of rj = nb/n y , the baryon to photon ratio, which is the only free parameter now left in the standard 
model. The information provided by Li is especially interesting in view of the recently conflicting observations of deuterium 
in high redshift absorption systems (Songaila et al. 1994, Carswell et al. 1996, Wampler et al. 1996, Rugers & Hogan 1996, 
Tytler, Fan & Buries 1996, Buries & Tytler 1996) . 

In the recent years the existence of a true plateau in the Li abundances of the warm halo dwarfs has been debated. 
Deliyannis, Pinsonneault & Duncan (1993) claimed the existence of a rather small dispersion of 10%, at 1 a level, in Li 
abundances on the Spite plateau. Thorburn (1994) found an even larger dispersion of Li abundance on the plateau and found 
a trend both with T e g and [Fe/H]. Norris, Ryan & Stringfellow (1994) and Ryan et al. (1996) also support the existence 
of a trend of Li abundance with both T e ff and [Fe/H] . The existence of dispersion and trends on the plateau could argue in 
favour of depletion and/or Galactic enrichment in halo dwarfs. In particular, stellar models have been proposed where Li is 
considerably depleted by rotational mixing, diffusion or stellar winds (Vauclair & Charbonnel 1995). If this is the case then 
the abundance of Li in warm Pop II dwarfs does not reflect the primordial Li abundance. 

On the other hand the existence of the Spite plateau has been reasserted by Molaro, Primas & Bonifacio (1995a) from 
the analysis of a subsample of stars with uniform T e s derived from Balmer line profiles by Fuhrmann, Axer & Gehren (1994). 
Their result was that there is no dispersion on the plateau above what is expected by observational errors and that there is 
no slope in either T e g or [Fe/H] . Spite et al. (1996) tackled the problem in a similar fashion using three samples of stars for 
which T e g was determined from either excitation equilibria, Balmer line profiles, or (b — y)o, and concluded that there is no 
evidence of dispersion or trends. So it appears that the existence of the Spite plateau is real when small samples of stars with 
homogeneous T e s are examined, but it is blurred when the whole set of stars is considered. 

In this paper we use a new T e s scale based on the semi-direct Infrared Flux Method (IRFM, Blackwell et al. 1990) 
applied to a large sample of stars by Alonso et al. (1996a), to investigate the distribution of Li abundance with [Fe/H] and 
T e ff for a statistically significant fraction of the stars with Li measured. Preliminary results have been already presented by 
Bonifacio & Molaro (1996). 
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Li is present mainly in the form of the singly ionized ion in the atmospheres of cool dwarfs and a large ionization correction is 
required when deriving the Li abundance from the Li I 6707 A resonance doublet. This correction is strongly dependent on 
the adopted T e g making the determination of T e g crucial for the discussion of possible trends and/or dispersion on the Spite 
plateau. Other stellar parameters such as surface gravity, metallicity or microturbulence are much less important for the Li 
abundance. In this paper we consider a new sample of stars with Li observations available and with accurate and internally 
coherent T e fi's. 

Effective temperatures are usually determined either from the continuum spectrum (colours, flux distribution) or from the 
line spectrum (Balmer line profiles, excitation equilibria). These methods are calibrated on the few stars for which the direct 
measure of T ef f is obtained from the bolometric fluxes and angular diameter measurements. The Sun is the only star cooler 
than F5V for which the angular diameter is known determining an intrinsic uncertainty in the scale of effective temperatures 
in the low main sequence. In particular, for Pop II stars there are no stars with measured angular diameter which can be 
used as primary calibrators; therefore, one is forced to rely on theoretical models to derive the metallicity dependence of his 
favourite T e g indicator. 

In the impossibility of obtaining a direct measure of T e g , the semi-direct Infrared Flux Method (IRFM, Blackwell et 
al. 1990) is the best alternative to derive the temperature of F and G stars. The IRFM relies only weakly on theoretical models 
and its main uncertainty is related to the absolute flux calibrations of the infrared magnitudes. It has been extensively applied 
to different samples of Pop I stars (Saxner & Hammarback 1985, Blackwell et al. 1990). The works of Magain (1987) and 
Arribas & Martinez- Roger (1989) applied the IRFM to limited samples of metal poor stars. Alonso et al. (1996a) have recently 
determined IRFM temperatures for an enlarged sample of 475 dwarfs of spectral types from F5 to K0, with a wide range of 
metallicities. The effective temperatures are reddening corrected. Alonso et al. (1996a) derived three effective temperatures 
for the J,H and K fluxes showing a good internal consistency. They provided a final T c g as the mean of the three, weighted 
with the inverse of their errors. The zero point of the IRFM T c g by Alonso et al. (1996a) differ by 112, and -56 K from 
the zero point of the IRFM T c « derived by Magain (1987), Saxner & Hammarback (1985) and Bell & Gustafsson (1989) 
respectively, which may be explained by several improvements in the Alonso et al. (1996a) analysis. The mean accuracy of 
the T e ff estimated by Alonso et al. (1996a) is about 1.5 %, which includes both systematic and random errors. 

Out of the sample of Alonso et al. (1996a), 64 stars have been observed for Li, and they form our sample for the study of 
the Spite plateau. The star names together with the stellar parameters and Li equivalent widths are reported in Table 1. The 
Li EWs for the 64 stars have been taken from the literature. For multiple measurements we adopted the weighted average. 
The errors in the EW have been taken from the original papers, when available, or from Ryan et al. (1996), which estimated 
the errors according to the Cayrel (1988) prescriptions. Following Ryan et al. (1996) we kept a 1 mA as the highest precision 
claimed in the observations, with the only exception of the few stars observed by Deliyannis, Boesgaard & King (1995) with 
the Keck telescope, at very high S/N. Multiple observations are in general in excellent agreement and an accurate statistical 
analysis of the errors for multiple observations of the same stars has been performed by Ryan et al. (1996) confirming that 
in the majority of cases the scatter between multiple observations is consistent with the stated errors. The measurement of 
EWs may be occasionally affected by non gaussian noise, such as that arising from cosmic rays hits or scattered light in the 
spectrograph. This event is rare but has to be considered in case of claims of intrinsic dispersion of isolated stars for which 
only single observations are available. 

The frequency of binarity among halo dwarfs is not well established but it may be between 20-30% up to 50% as it is 
for Pop I stars. An investigation on the impact of binarity was done by Molaro (1991) showing that known binaries were 
not significantly different from other stars. However, the influence of binaries may be important in assessing the presence of 
a small dispersion. In Molaro et al. (1995a) the A(Li) in three stars was consistent with the mean value only considering 
the full error. One of these, HD 116064, has been later found to be a spectroscopic binary by Spite et al. (1996) and the 
compensation for the veiling, which is of 20 % in the EW, move the A(Li) closer to the mean value. This is a clear example 
of how binarity may affect the dispersion analysis on the plateau. Another example is G020-008 found to be double by Fouts 
(1987), which, in the small sample of stars considered by Spite et al. (1996), contributes to increase the dispersion. 

Several binaries are present in our sample. HD 219617, BD+20 3603, HD 3567, HD 16031, BD+01 2341p, HD 84937, 
HD 132475, G206-034, HD 188510, BD+38 4955 have been identified as binaries or suspected radial velocity variables by 
Carney (1983). By means of spectroscopic analysis over a four-year temporal baseline, Stryker et al. (1985) classified as 
certain binaries: HD 84937, HD 94028, HD 201891 BD+26 2606, G090-025, HD 188510; as possible binaries: BD+28 2137, 
BD+34 2476, BD+13 3683, BD+ 26 3578; near significance criterion: BD+01 2341p, BD+29 2091, BD+38 4955, BD+21 607, 
HD 108177, G090-003. Lu et al. (1987) by means of the speckle technique identified: BD+17 4708. G020-008, G020-024, HD 
16031 have been found to be binaries by Fouts (1987). Thus in total there are 25 known or suspected binaries out of 64 stars, 
which gives a binarity fraction of at least 39% in our sample. 
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Figure 1. a) cl, (b — y) diagram for our sample of Pop II stars; b) cO, (b — y)o diagram. The solid line represents the locus of points with 
log g = 3.5 for [Fe/H]=-1.5. The dashed line is the same but for [Fe/H]=-3.0. The region above these lines is populated by subgiants and 
giants. The crosses are the stars with [Fe/H] > —1.5 while the X symbols are stars with [Fe/H] < —1.5. 

3 MODEL ATMOSPHERES AND Li ABUNDANCES 

The different atmospheric models adopted by different authors to derive lithium abundances are known to be generally 
consistent with one another, but in some cases they may be responsible for systematic differences in lithium abundances. 
As shown by Molaro et al. (1995a) and Molaro, Bonifacio & Primas (1995b), the Kurucz (1993) models which include 
overshooting are hotter than the Bell & Gustafsson models in the Li doublet forming region. As a consequence, the Thorburn 
(1994) Li abundances, derived by using Kurucz models with overshooting, are systematically « 0.1 dex higher than those of 
other researchers. On the other hand the Kurucz (1993) models with the overshooting switched off are in good agreement 
with the others. Castelli, Gratton & Kurucz (1996) have shown that for all F and G stars, except the Sun, the Kurucz models 
computed without overshooting are more consistent with the observations than models with overshooting. Consistently with 
these findings, in this investigation we use model atmospheres computed with the ATLAS9 code (Kurucz 1993) with the 
overshooting option switched off. 

For the opacities we used the ODFs with enhanced a elements, which provide a more realistic chemical composition for 
Pop II stars than the solar scaled ODFs. The microturbulent velocity was assumed to be 1 kms -1 . However, the precise value 
of the microturbulent velocity is relatively unimportant for the Li abundance in halo stars, since 0.5 kms _1 change in the 
microturbulent velocity induce a change of 0.005 dex in Li. 

We checked our models with those used in Ryan et al. (1996) by comparing our curve of growth for a model with 
T e fi = 6500 , logg= 4.0, [Fe/H] -2.0 and microturbulence of 1 kms - with the corresponding curve of growth, based on Bell- 
Gustafsson models, published by Ryan et al. (1996). The difference in Li abundance is zero for an equivalent width of 14 
mA (A(Li)w 2.1) and 0.008 dex for an equivalent width of 24 mA (A(Li)ss 2.4), with our models giving abundances higher 
than the Bell-Gustafsson ones. 

Surface gravities have been redetermined for each star to identify possible subgiants. As already mentioned, the Li doublet 
equivalent width depends little on surface gravity with a 0.02 dex change in Li abundance for a 0.7 dex change in logg, but 
the identification of subgiants in the sample is important because their Li abundance is expected to be lower due to dilution 
(Deliyannis et al. 1990). As a gravity indicator we used the cO index of Stromgren photometry, which measures the Balmer 
discontinuity. We derived surface gravities from cO in the same way as described in Molaro et al. (1996). For all our stars, 
except G201-005, we took the photometry from the General Catalogue of Photometric Data (Mermillod, Hauck & Mermillod 
1996), this was dereddened using the Schuster & Nissen (1989) intrinsic colour calibration. For BD +71 31, for which there is 
no j3 available, we assumed a zero reddening. For some stars we could not find any solution; this happened if the observed cO 
value was smaller than any of the theoretical values for the given temperature and metallicity. This was the case for G246-038, 
HD25329, G090-025 and G190-015. For these stars, as well as for G201-005, we have taken the surface gravity given by Alonso 
et al. (1996a). From our model grid, described in Molaro et al. (1996), we computed fluxes and both Johnson and Stromgren 
photometry. In Fig. 1 the cl, {b — y) and cO, (b — y)o diagrams for our sample of Pop II stars are shown. The solid line represents 
the locus of points with log g = 3.5 for [Fe/H]=-1.5, from our synthetic Stromgren photometry, and the dashed line is the 
same but for [Fe/H]=-3.0. The region above these lines is populated by subgiants and giants. In our sample of 41 plateau 
stars only 3 have a surface gravity less than or equal to 3.6, namely G090-003, G126-062 and BD+71 31, with logg 3.55, 3.58 
and 3.58, respectively. 
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Figure 2. The A(Li)- [Fc/H] diagram for our sample of stars. The circles are stars with T e g > 5700 K, while the squares are those with 
^eff < 5700 K. Subgiants are shown with crossed symbols. Three upper limits are shown as vertical downward arrows. 

A model atmosphere was computed for each star with the appropriate gravity and metallicity. We then computed a 
synthetic profile of the Li doublet with the SYNTHE code (Kurucz 1993) and obtained the equivalent width of the doublet 
by integrating this profile. The process was iterated, changing the Li abundance, until the computed equivalent width matched 
the observed one. This approach is different from that used by Thorburn (1994), Molaro et al. (1995a) and Spite et al. (1996), 
who treated the doublet as a single line. The differences are very small, because for our stars the stronger Li line is never 
saturated, however the explicit computation of the doublet is physically more realistic. 

The present analysis is carried out strictly under the LTE assumption. NLTE effects have been studied by Carlsson et 
al. (1994) and Pavlenko & Magazzu (1996) and were shown to be minimal for halo dwarfs. The corrections provided by 
Carlsson et al. (1994) have been applied to the LTE abundances and considered in the discussion of the results. 

The errors reported in Table 1 are those following from the error in EW and in T c ff . The two errors are summed under 
quadrature as it is generally done. The errors given do not consider a possible error originated by uncertaintes in the gravity 
and in the microturbulence which are of the order of 0.01 dex each. These are also independent sources of errors and should 
be added to the above quoted errors. We found that ignoring their contribution to the total error does not affect our results. 

4 RESULTS 

4.1 Li - T of f ,[Fe/H] correlations 

The derived Li abundances are given in Table 1 (where A(Li) = log(A?Li/./Vff ) + 12). The corresponding A(Li) - T e s and A(Li) 
- [Fe/H] diagrams are given in figures 2 and 3. We considered plateau stars those with T e tt > 5700 and with [Fe/H] < —1.5; 
this extracts a sample of 41 stars. 

To investigate the existence of trends in the A(Li) - T e g and A(Li) - [Fe/H] planes, we first performed univariate fits by 
means of four different estimators: (i) the Bivariate Correlated Errors and intrinsic Scatter (BCES) of Akritas & Bershady 
(1996); (ii) BCES simulations bootstrap based on 10000 samples; (iii) least squares with errors in the independent variable only 
(as implemented in the routine f itxy of Press et al. (1992)); (iv) least squares with errors in both variables (as implemented 
in the routine f itexy of Press et al. (1992)). In addition, the presence of a slope was searched by means of a rank correlation 
analysis using Kendall's r. As a second step we performed bivariate fits using both T e ff and [Fe/H] as independent variables. 
The need for such bivariate fits follow from the known trend of T e fi with [Fe/H] , which is an observational bias, due to the fact 
that the more metal-poor stars are, on average, the more distant and therefore also the more luminous and hotter members 
of the class. So one could suspect that if there is, indeed, a trend of A(Li) with T e g and [Fe/H] , the trend of T c g with 
[Fe/H] could be exactly tuned so as to suppress the trends when performing univariate fits. For the bivariate fits we used an 
ordinary least squares method taking into account only errors in the independent variable. 

The results of the statistical analysis are reported in Tables 2-4. Inspection of Table 2 shows that a small slope of the 
order 0.02(± 0.01)/100 K is marginally detected using BCES and fitxy, but not using f itexy. We used Kendall's r to check 
the detection of the slope and we found a r = 0.2112 with a two-sided significance level of 0.05, consistent with a positive 
detection of the slope. Note that the slope found is about one half of those found by Thorburn (1994) of 0.034(± 0.006)/100 K 
and Ryan et al. (1996) of 0.0408(± 0.0052)/ 100K. One may also note that the goodness-of-fit implied by fitxy is quite small, 
i.e. of the order of 0.035. This value means that the fit may be believable if the errors are nonnormal or if they have been 
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Table 1. Li data 



BD HD G# T cB aT cB logg [Fe/H] EW oEW A(Li) A(Li)i A(Li) 2 A(Li) 3 a T a EW a* u 



BD 


— 09 


122 b 


HD 3567 


0270-023 


5858 


63 


3.70 


-1.34 


45.0 


6 





2 


28 


2, 


.30 


2. 


32 


2, 


.33 


0.05 


0.08 


0.10 


BD 


4- 71 


31 






6026 


94 


3.58 


-1.91 


31.0 


2 


3 


2 


19 


2 


20 


2 


.21 


2. 


22 


0.07 


0.04 


0.08 


BD 


-1- 29 


366 




0074-005 


5685 


85 


4.48 


-1.05 


14.0 


2 


9 


4 


55 


4 


57 


1 


.64 


4 


66 


0.07 


0.11 


0.13 


BD 


— 13 


482 b 


HD 1 6031 




6114 


72 


4.00 


-2.09 


28.0 


2 


2 


2 


20 


2 


21 


2 


.21 


2. 


23 


0.05 


0.04 


0.06 


BD 


+ 09 


352 




O076-021 


5931 


85 


4.03 


-2.15 


34.0 


I 


.5 


2 


16 


2 


.18 


2 


.21 


2 


.22 


0.06 


0.03 


0.07 










0004-037 


6337 


92 


4.13 


-3.31 


20.0 


2 


2 


2 


45 


2 


16 


2 


.15 


2 


16 


0.06 


0.05 


0.08 


BD 


+ 25 


495 


HD 1 9445 


O037-026 


6050 


109 


4.45 


-2.15 


34.0 





5 


2 


26 


2 


.27 


2 


28 


2 


30 


0.08 


0.01 


0.08 


BD 


+ 66 


268 






5279 


81 


4.50 


-2.94 


10.0 




g 


4 


09 


4 


44 


1 


74 


4 


76 


0.08 


0.65 


0.66 


BD 


— 04 


680 


HD 24289 




5866 


87 


3.73 


-2.07 


47.0 


2 


4 


2. 


29 


2 


31 


2 


35 


2 


36 


0.07 


0.03 


0.08 








HD 25329 




4842 


97 


4.50 


-1.64 


< 3.0 






o 


00 




















BD 


4- 21 


607 b 


HD 284248 


O008-01 6 


6034 


74 


4.05 


-1.61 


25.0 


2 


2 


2 


44 


2 


12 


2 


42 


2 


.14 


0.05 


0.05 


0.07 


BD 


+ 03 


740 




0084-029 


6110 


82 


3.73 


-2.01 


21.0 




.0 


2 


06 


2 


07 


2 


.07 


2 


09 


0.06 


0.02 


0.06 


BD 


4- 37 


1458 




0098-058 


5510 


56 


3.80 


-1.86 


14.0 




4 


4 


41 


4 


43 


1 


.65 


1 


67 


0.05 


0.05 


0.07 










0088-01 

VJUOO U 1U 


5856 


74 


3.85 


-2.20 


30.0 


2 


.9 


2 


05 


2 


07 


2 


.12 


2 


13 


0.06 


0.05 


0.08 






b 




0090-003 


5786 


132 


3.55 


-1.99 


44.0 




9 


2 


.19 


2 


21 


2 


.27 


2 


29 


0.10 


0.03 


0.10 


BD 


4- 24 


1676 




VjjUOO uoz 


6208 


67 


3.78 


-2.59 


27.0 


l 


g 


2 


23 


2 


24 


2 


25 


2 


26 


0.05 


0.04 


0.06 


BD 


4- 31 


1684 b 


HD 64090 




5441 


68 


4.50 


-1.82 


12.0 


I 


.0 


4 


30 


4 


32 


1 


.61 


4 


.63 


0.07 


0.04 


0.08 


BD 


— 15 


2546 


HD 74000 




6224 


96 


4.08 


-2.00 


25.0 


3 


4 


2 


22 


2 


23 


2 


.23 


2 


.24 


0.07 


0.07 


0.10 


BD 


— 15 


2656 


HD 76932 




5727 


77 


3.78 


-1.05 


28.7 




.0 


4 


.94 


4 


96 


2. 


01 


2. 


03 


0.06 


0.02 


0.06 










vjj _l iu-uiy 


5950 


72 


5.00 


-1.86 


38.0 


3 


4 


2 


26 


2 


28 


2 


29 


2 


30 


0.05 


0.05 


0.07 


BD 


+ 09 


2190 




O041 -041 


6333 


89 


3.83 


-2.89 


19.0 




4 


2 


44 


2 


15 


2 


.14 


2 


16 


0.06 


0.04 


0.07 


BD 


4- 01 


2341p b HD 83769C 1 


0048-029 


6313 


80 


3.90 


-2.81 


22.0 


2 


3 


2 


20 


2 


.21 


2 


.21 


2 


.22 


0.06 


0.05 


0.08 


BD 


4- 14 


2151 b 


HD 84937 


0043-003 


6330 


83 


3.90 


-2.49 


24.0 




2 


2. 


26 


2 


.27 


2 


.26 


2, 


.28 


0.06 


0.03 


0.07 


BD 


4- 29 


2091 b 




Ol 1 9-032 


5815 


119 


4.80 


-1.86 


45.0 


g 


9 


2 


24 


2. 


26 


2 


30 


2. 


32 


0.09 


0.09 


0.12 


BD 


4- 21 


2247 b 


HD 94028 


0058-025 


6001 


131 


4.15 


-1.50 


35.0 




.6 


2. 


25 


2. 


26 


2 


.27 


2 


28 


0.10 


0.02 


0.10 


BD 


+ 02 


2406 


HD 97916 




6393 


155 


3.98 


-1.22 


< 3.0 






4 


.28 


4 


29 


4 


.32 


1 


33 


0.10 














O01 1 -044 

V. 1 V ' 1 _1_ \J ' L 'J. 


5943 


108 


4.08 


-1.80 


29.0 


3 


4 


2 


40 


2 


12 


2 


.13 


2 


.14 


0.08 


0.06 


0.10 


BD 


+ 02 


2538 b 


HD 108177 


O01 3-035 


6097 


94 


4.03 


-1.98 


32.0 




3 


2 


.26 


2 


.27 


2 


.27 


2. 


.29 


0.07 


0.02 


0.07 










0059-024 


6108 


83 


4.28 


-2.69 


36.0 


2 




2 


31 


2 


.32 


2 


.34 


2 


35 


0.06 


0.03 


0.07 


BD 


+ 28 


2137 b 




0059-027 


6107 


127 


3.85 


-2.34 


27.0 


2 


2 


2 


17 


2 


18 


2 


.19 


2 


.20 


0.09 


0.05 


0.10 










0064-01 2 


6468 


87 


4.10 


-3.35 


27.0 


x 


.8 


2 


39 


2 


40 


2 


.36 


2 


37 


0.06 


0.03 


0.07 


BD 


4- 34 


2476 b 




ri 1 fir: non 


6282 


72 


3.88 


-2.41 


23.0 


I 


.5 


2 


21 


2 


22 


2 


22 


2 


23 


0.05 


0.04 


0.06 










VJUUt: UO 1 


6432 


70 


4.15 


-2.51 


15.0 


I 


.0 


2 


40 


2 


.11 


2 
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Table 1. Li data (continued) 



BD HD 


G# 


T c ff 


°~T eS 


logg 


[Fe/H] 


EW 


uEW 


A(Li) 


A(Li)i 


A(Li) 2 


A(Li) 3 


erf 


&EW 


Ctoi 


BD + 02 4651 
BD + 59 2723 


G029-023 
G217-008 


6102 
6134 


90 
70 


3.80 
4.25 


-2.02 
-1.91 


27.0 
29.0 


2.2 
2.1 


2.18 
2.23 


2.19 
2.24 


2.19 
2.24 


2.21 
2.25 


0.07 
0.05 


0.05 
0.04 


0.09 
0.06 



A(Li)i is the NLTE corrected value according to Carlsson et al. (1995) 

A(Li)2 is the value corrected for the depletion predicted by the standard Li isochrones of Deliyannis et al. (1990) 
A (Li) 3 is the value corected both for NLTE and standard depletion 

* the sum under qaudrature has been made taking into account three significant figures and then rounded off. 
b binaries 

1 this star is LP 608-62=Ross 889, it is indicated as HD 83769C by Olsen (1993) 



Figure 3. The A(Li) - T c g diagram for our sample of stars. The circles are stars with [Fc/H] < —1.5 , while the squares are those with 
[Fc/H] > —1.5. Subgiants are shown with crossed symbols. Three upper limits are shown with vertical arrows pointing downwards. 

slightly underestimated (Press et al. 1992). To check this interpretation we performed the same fits, but using the slightly 
larger errors which you obtain by summing linearly, rather than under quadrature, the errors due to T e g and EW. While 
the value of the fitting parameters are little sensitive to this assumption, the goodness-of-fit boosts up to values of the order 
of 0.7. This comment applies also to the results of Table 3, obtained from f itxy, and to those of Table 4. 

Table 3 shows the results of the investigations of trends of A (Li) with [Fe/H]. Irrespective of the method used the slope 
is not detected, Kendall's r is -0.857 with a two-sided significance level of 0.43, i.e. no slope. Note that in almost all cases 
the slope with [Fe/H] is negative, contrary to the findings of Thorburn (1994) and Ryan et al. (1996). A negative slope is 
physically plausible since it could be present if astration is important. An increase of Li with the metallicity is expected when 
Li synthesis by GCR begins to become significant. Using Be as an indirect indicator, from the Be-metallicity correlation 
found by Molaro et al. (1996), the contribution to Li by energetic galactic cosmic rays is of only 0.03 dex at [Fe/H]=-1.5 and 
of 0.1 dex at [Fe/H]=-1.0, which is below our detectability threshold. 

Table 4 displays the results for the bivariate fits, which show that the values of the slopes found by performing bivariate 
fits are not significantly different from the values found performing univariate fits. 

According to the Deliyannis et al. (1990) standard ZAMS Li isochrones, a very little depletion is indeed expected in the 
T e ff range between 5700 K and 6400 K. In Fig. 4 one can see the Li isochrones of Deliyannis et al. (1990) corresponding 
to [Fe/H] =— 2.3 and [Fe/H] = — 1.3, and an age of 16 Gyr, which have been shifted upwards by 0.1 dex to better match the 
observations. Using these isochrones we corrected the observed Li abundances so that each star was assigned the A(Li) it 
should have if its temperature were 6400 K and its metallicity -2.3. In Fig. 5 we show the corrected data and no trend is 
evident. In fact the formal statistical analysis of this corrected data shows no detectable trend irrespective of the estimators 
used. Kendall's r is 0.0804 with a two-sided significance level of 0.458, consistent with a non-detection of the slope. Thus 
the small slope with T e s , which is marginally detected using our sample, is consistent with a constant Li abundance on the 
plateau with minimal depletion predicted by the standard Li isochrones of Deliyiannis et al. (1990). Also for the case of 
bivariate fits the small slope, marginally detected using the uncorrected Li, disappears after accounting for the Li depletion 
predicted by standard isochrones. 

NLTE correction are even smaller, but again they act in the sense of decreasing the slope. Thus when both corrections 
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Table 2. Univariate fits in the A(Li)-T e g plane. 



Method 



A(Li) = 0.766(±0.635) + 0.0235(±0.0104) x (T e// /100) 
A (Li) = 0.774(±0.644) + 0.0233(±0.0106) x (T e// /100) 
A (Li) = 1.346(±0.353) + 0.0140(±0.0058) x (T e// /100) 
A (Li) = 1.153(±4.056) + 0.0171(±0.0511) x (T e// /100) 



BCES 

BCES bootstrap 



fitxy x 



56.59 ; P = 0.034 



f itexy x = 0.4291 ; P = 1.000 



after the theoretical correction for the standard depiction 



A(Li) = 1.63 (±0.637) 
A (Li) = 1.64 (±0.646) 
A (Li) = 2.009(±0.353) 
A (Li) = 1.827(±4.053) ■ 



0.0097(±0.0104) x (T e// /100) 
0.0095(±0.0106) x (T e// /100) 
0.0035(±0.0058) x (T e// /100) 
0.0064(±0.0444) x (T e// /100) 



BCES 

BCES bootstrap 

fitxy x 2 = 56.47 ; P = 0.035 

f itexy x 2 = 0.4250 ; P = 1.000 



after the theoretical correction for the NLTE effect 



A(Li) = 0.858(±0.635) + 0.0222(±0.0104) x (T e// /100) 
A (Li) = 0.866(±0.644) + 0.0220(±0.0106) x (T e// /100) 
A (Li) = 1.421(±0.353) + 0.0130(±0.0058) x (T e// /100) 
A (Li) = 1.230(±4.056) + 0.0160(±0.0503) x (T e// /100) 



BCES 

BCES bootstrap 



fitxy x 
f itexy x 2 



56.59 ; 
0.4291 



0.034 
= 1.000 



after the theoretical corrections for the standard depletion and NLTE 



A(Li) = 1.72 (±0.637) + 0.0084(±0.0105) x (T e// /100) 
A(Li) = 1.73 (±0.646) + 0.0083(±0.0106) x (T e// /100) 
A (Li) = 2.084(±0.353) + 0.0025(±0.0058) x (T e// /100) 
A (Li) = 1.902(±4.053) + 0.0054(±0.0441) x (T e// /100) 



BCES 

BCES bootstrap 
fitxy x 2 = 56.47 ; P = 
f itexy x 2 = 0.4250 ; P 



0.035 
= 1.000 



Figure 4. Zoom of the A(Li) — T e fj diagram for our sample of stars. The solid line is the ZAMS isochronc of Dcliyannis ct al. (1990) 
corresponding to [Fe/H] = —1.3 and an age of 16 Gyrs, the dashed line is the same but for [Fe/H] = —2.3. Both isochroncs have been 
shifted upwards by 0.1 dex to better match the observations. 

are applied no slope is detected. 

As far as the slope is concerned, a possible explanation of the different results may lie in the different T e ff adopted. 
Thorburn (1994) derived a common photometric temperature scale by adapting the temperature calibrations between T c g - 
(V-K); (b-y) of Carney (1983) and Carney et al. (1987) into a T e g -(B-V) scale, which is assumed to be valid for all metal 
poor stars. There are 24 stars in common with Thorburn (1994) and in Fig. 6a the temperature difference between Alonso 
et al. (1996a) and Thorburn (1994) shows a systematic trend with T e H , with the T e g of Alonso et al. (1996a) higher at 
the hot edge of the plateau and lower at the cool edge. This systematic difference in the T c s produces the slope in Li versus 
T e ff which is observable in Fig. 6b, where the differences between our A(Li) values and those of Thorburn (1994) are shown 
as a function of T e a . In addition to the trend one may notice an offset, which is due to the inclusion of overshooting in the 
ATLAS9 models used by Thorburn. All our stars, but one, are in common with the sample of Ryan et al. (1996) and in Fig. 
6c and 6d we compare our results with those of Ryan et al. (1996) in the same way we did for the Thorburn data. Ryan 
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Figure 5. Zoom of the plateau in the A(Li) — T c ff plane using the Li abundances corrected for standard depletion. The solid line is our 
estimate of the primordial Li, the two dashed lines are the la errors. 

Table 3. Univariate fits in the A(Li)-[Fe/H] plane. 



Method 



A(Li) 
A(Li) 
A(Li) 
A(Li) 


= 2.09(±0.09) - 0.05(±0.04) x [Fe/H] BCES 

= 2.10(±0.09) - 0.04(±0.04) x [Fe/H] BCES bootstrap 

= 2.15(±0.06) - 0.02(±0.03) X [Fe/H] f itxy X 2 = 61.49 ; P 

= 2.11(±0.15) - 0.04(±0.06) X [Fe/H] f itexy X 2 = 8.50 ; P = 


= 0.012 
1.000 


after the theoretical correction for the standard depletion 


A(Li) 
A(Li) 
A(Li) 
A(Li) 


= 2.17(±0.08) - 0.02(±0.04) x [Fe/H] BCES 

= 2.18(±0.08) - 0.02(±0.03) x [Fe/H] BCES bootstrap 

= 2.22(±0.06) - 0.00(±0.03) X [Fe/H] f itxy X 2 = 56.84 ; P 

= 2.18(±0.15) - 0.02(±0.06) X [Fe/H] f itexy X 2 = 7.67 ; P = 


= 0.032 
1.000 


after the theoretical correction for the NLTE effect 


A(Li) 
A(Li) 
A(Li) 
A(Li) 


= 2.11(±0.09) - 0.04(±0.04) x [Fe/H] BCES 

= 2.12(±0.08) - 0.04(±0.04) x [Fe/H] BCES bootstrap 

= 2.17(±0.06) - 0.02(±0.02) X [Fe/H] f itxy X 2 = 60.90 ; P 

= 2.13(±0.15) - 0.03(±0.06) X [Fe/H] f itexy X 2 = 8.41 ; P = 


= 0.014 
1.000 


after the theoretical corrections for the standard depletion and NLTE 


A(Li) 
A(Li) 
A(Li) 
A(Li) 


= 2.19(±0.08) - 0.02(±0.04) x [Fe/H] BCES 

= 2.20(±0.08) - 0.02(±0.03) x [Fe/H] BCES bootstrap 

= 2.24(±0.06) + 0.00(±0.02) X [Fe/H] f itxy X 2 = 56.66 ; P 

= 2.20(±0.15) - 0.01(±0.06) X [Fe/H] f itexy X 2 = 7.64 ; P = 


= 0.033 
1.000 



Table 4. Bivariatc fits 











x 2 


P 


Li correction 


1.10(±0.43) + 0.034(±0.034) 


X 


[Fe/H] + 0.019(±0.008) 


x (T e// /100) 


60.21 


0.012 


none 


1.88(±0.43) + 0.018(±0.034) 


X 


[Fc/H] + 0.006(±0.008) 


x (T e// /100) 


61.03 


0.010 


standard depletion 


1.17(±0.43) + 0.034(±0.034) 


X 


[Fc/H] + 0.018(±0.008) 


X (Te///100) 


60.57 


0.011 


NLTE 


1.95(±0.43) + 0.018(±0.034) 


X 


[Fc/H] + 0.005(±0.008) 


x (T e// /100) 


61.14 


0.010 


standard depletion and NLTE 



et al. (1996) adopted the average of the {b-y)-T e s from Magain (1987), the (R-I)-T e ff from Buser & Kurucz (1992) and a 
(B-V)-T c £f worked out by interpolation from a grid obtained from the literature. Also in this case, as can be seen in Fig. 6d, a 
temperature scale difference, which is responsible for the differences in the slope A(Li)-T e £f found in the two analyses, is clear. 
Although the T e g scales used by Thorburn (1994) and Ryan et al. (1996) may be sufficient for a general abundance analysis, 
they may be not adequate to discuss dispersion and/or trends at the levels claimed, since the indices used (b — y, B — V and 
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R — I) are all metallicity dependent and the first two depend also on gravity. To illustrate the effect of metallicity we divided 
the 40 stars in common with Ryan et al. (1996) in four metallicity bins (-2.0 < [Fc/H] < -1.5, -2.5 < [Fe/H] < -2.0 , 
—3.0 < [Fe/H] < —2.5 and [Fe/H] < —3.0), the bins have 16, 13, 8 and 3 stars respectively. For each bin we compute the 
mean of Taioiiso TRy an , the results (rounded to the nearest degree) are: +10 K, +50 K, +96 K and +132 K. 

The same analysis with the 16 stars we have in common with Molaro et al. (1995a) shows a considerable scatter, but no 
obvious trend. They used the T c g obtained by Fuhrmann et al. (1994) by fits of the Balmer lines, i.e with a totally different 
approach. Balmer line temperatures have the considerable advantage that they are reddening independent. A trend with 
temperature cannot be excluded, but the number is too small to draw any firm conclusion. Ryan et al. (1996), reanalysed the 
Molaro et al. (1995a) sample, after excluding the two subgiants and outliers, and found a trend of « 0.04/100 K. However, 
no trends are found using the whole original sample. 

We thus conclude that the presence or absence of trends of lithium abundance with T e g appears to be strongly dependent 
on the temperature scale adopted. 

4.2 Dispersion on the plateau? 

The straight average value of A(Li) for the selected stars is A(Li) = 2.19 with a dispersion around the mean value of ±0.094. 
This dispersion is compatible with the root mean square of the estimated observational errors (0.083 dex) showing no clear 
evidence for intrinsic dispersion. In fact 29 stars (« 71%) are within 1 a of the mean A (Li), which is compatible with a normal 
distribution, such as would be expected for observational errors. The dispersion around the mean value further decreases 
down to 0.088 if the Li abundances corrected for the theoretical depletion and the NLTE effects are considered. In Fig. 5 
the corrected Li abundances on the plateau are zoomed illustrating that the data are consistent with a unique Li abundance 
within the errors. The absence of intrinsic scatter is remarkable in the light of the high percentage of binaries or suspected 
binaries in our sample. 

Thorburn (1994) analysed the presence of scatter in the equivalent width - T e g plane, claiming that the null hypothesis 
of no intrinsic dispersion is rejected at the 6a confidence level. In Fig. 8a we show such a diagram for our plateau sample, 
together with the theoretical equivalent width - T e g locus obtained from our curves of growth corresponding to three Li 
abundances (2.20 ± O.ldex) for a metallicity of [Fe/H]=-2.0. Note that this locus is not a straight line as assumed both by 
Thorburn (1994) and Deliyannis et al. (1993), but is better approximated by an exponential, as can be seen by the linear 
relation displayed in the \og(EW)-T e g diagram shown in Fig. 8b. The slope of the theoretical lines is w —0.0623 log(mA)/ 
100 K. Performing a simple least squares fit with a straight line of this slope (i.e. fitting only one free parameter) we find the 
dotted line shown in Fig 8b. The fit has ax 2 — 0.393, corresponding to a goodness-of-fit of practically 1.00, and an R.M.S 
= 0.066 log(mA). The fitting line almost coincides with the theoretical line corresponding to A(Li) = 2.20. This shows that, 
even neglecting the metallicity effects on the curve of growth, the data are consistent with the hypothesis of a constant Li, 
around 2.20, and a scatter due to observational errors only. 

The disagreement with Thorburn (1994) on this point can be explained by the fact that she fits an exponential with a 
straight line. She attributes the very low goodness of fit (« 10~ 10 ) to either an intrinsic scatter in the Li abundances or to an 
underestimation of the errors by about 50 %, while this can be explained by the use of an incorrect fitting function. In addition, 
as pointed out by Spite et al. (1996), the temperature errors of Thorburn (1994) appear to be somewhat underestimated. 

Deliyannis et al. (1993) discussed the dispersion on the Spite plateau in the W\ — (b — y)o plane which is strictly 
observational and avoided the problems implied in the determination of T e g . Note, however, that like Thorburn (1994) they 
also assumed a linear relation between W\ and (6 — y) to look for dispersion. Moreover, there is not a one-to-one mapping 
from (6 — y)o to T e g since this index depends both on gravity and metallicity. A sample of halo stars with different gravities 
and metallicities will show a scatter in the W\ — (b — y)o plane even if there is no scatter in the Li abundances, simply because 
when comparing stars of the same (6 — y)o we compare stars with different T cH . 

Ryan et al. (1996) identified a triplet of stars (G064-012, G064-037, and CD -33 1173) with similar metallicity and 
temperature, but significantly different Li abundances. Only the first two belong to our sample and in fact show one of the 
maximum differences in the A(Li) values in the whole sample, namely 2.39+0.07 and 2.10 ±0.06 , which agree at the 2.23 
a level. For comparison the values of Ryan et al. (1996) are 2.29±0.05 and 2.01 ±0.04. Our errors associated with the 
T e ff are in general larger than those of Ryan et al. (1996). The temperature errors of Ryan et al. (1996) are rather small 
(a(B — V) = 50K, a(b — y) — 35 K and a(R — I) = 90 K). If we compare them with the standard deviations of the colour- 
temperature calibrations derived by Alonso Arribas & Martinez-Roger (1996b) (a(B — V) = 130K, a(b — y) = 110 K and 
a(R — I) = 135 K) we see that the latter are larger by almost a factor of three. The calibrations of Alonso et al. (1996b) are 
based on over 400 stars of which roughly 50% are Pop II, while those of Magain (1987), on which Ryan et al. (1996) partially 
rely, are based on a total of 44 stars out of which only 11 belong to Pop II. This comparison gives rise to the suspicion that 
the temperature errors of Ryan et al. (1996) are underestimated. 

The absence of intrinsic dispersion found by us is in agreement with the results of Molaro et al. (1995a) and Spite et 
al. (1996). The latter authors have considered three small samples of stars with accurate T e g derived either from excitation 
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Figure 6. Temperature and A(Li) differences between the present paper and those of Thorburn (1994)(panels a) and b)), Ryan ct 
al. (1996) (panels c) and d)) and Molaro ct al. (1995a)(panels e) and f)), for the stars in common. 
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Figure 7. The temperatures of Thorburn (1994), panel a), Ryan et al. (1996), panel b), and Fuhrmann ct al. (1994), panel c), are 
plotted as a function of the temperatures of Alonso et al. (1996). The bisector is shown as a solid line. 



Figure 8. Equivalent widths as a function of effective temperature. The solid lines are the theoretical loci for A(Li)=2.20 and 
[Fe/H] = —2.00; the dashed lines are the same but for A(Li)= 2.20 ± 0.1. The dotted line in panel b) is the fit to the data of a 
line of slope —0.0623 log(mA)/100 K, which is the slope of the theoretical lines in the same panel. 

temperature or from the dereddened colour (b-y)o or from Ha wings. They found that the dispersion of the Li abundances in 
these small subsamples is significantly reduced, with an R.M.S. of the order of 0.06-0.08 dex, and they argue that it is fully 
explained by the errors in the T e ff and EWs. It is interesting that the presence of the binary G020-008, with a A(Li) lower 
than the average, is found to contribute considerably to the scatter. 

The absence of dispersion on the plateau is remarkable considering the number of processes that might produce it, 
such as localized Galactic enrichment, chromospheric Li enrichment, different amounts of accretion or astration, binarity, 
anomalous reddening, mistaken subgiants, pre main-sequence depletion and others. This implies that such possibilities are 
either intrinsically or statistically irrelevant. 



5 DISCUSSION 
5.1 Li depletion? 

We have found that in our sample of 41 stars on the plateau with the IRFM T e a derived by Alonso et al. (1996a) there is 
no evidence for intrinsic dispersion in excess of what can be expected from the measurement uncertainties. In addition, there 
is no evidence of trends with the metallicity, although our stars span more than two orders of magnitude in [Fe/H]. The tiny 
trend with T e fj , which has been detected, can be entirely explained by the standard Li isochrones of Deliyannis et al. (1990). 
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The view that the Li abundance on the Spite plateau is the primordial abundance requires that Li has not been destroyed, 
in appreciable amounts, in the atmospheres of halo dwarfs. Simple models of Li evolution, referred to as standard models, 
predict little or no depletion of Li on the MS, thus supporting the primordial nature of Li in halo dwarfs. However, it has 
been suggested that three different mechanisms are able to deplete Li in halo dwarfs in significant amounts. They are diffusion 
(Michaud, Fontaine & Beaudet 1984), rotational mixing (Pinsonneault, Deliyannis & Demarque 1992) and stellar winds 
(Vauclair & Charbonnel 1995). More elaborate scenarios with some combination of the three have been also considered. To 
reproduce the observations such depletion models must be able to deplete Li uniformly in a way which is independent from 
the metallicity and from the mass (and therefore T e g ) of the star. The detection of sizeable dispersion on the plateau would 
be a strong evidence in favour of some Li depletion during the life of halo stars. 

Diffusion has been proposed to be responsible for abundance anomalies in Ap and Am stars which are somewhat hotter 
than halo dwarfs, have non convective atmospheres and the Ap are strongly magnetic. In halo stars there are no abundance 
anomalies which have been ascribed to diffusion. In fact Li is the only element for which diffusion has been claimed to be 
effective in halo dwarfs. First computations of the effects of diffusion on the Li abundances in halo stars have been performed 
by Michaud et al. (1984). Diffusion causes Li to sink below the photosphere and Michaud et al. (1984) found that the 
process is more efficient at the hot end of the plateau producing a pronounced downturn of the Li abundance for T c g > 
6000 K. The depth of the surface convection zone is sensitive to opacities and to the mixing length, so that some model 
dependence is present. However the downturn persists in the various assumptions (Chaboyer & Demarque 1994, Vauclair & 
Charbonnel 1995). From an initial A(Li)=2.5, at T e g = 6500 K the predicted Li abundance becomes A(Li) = 1.5, while our 
observed abundance is 0.72 dex higher. The total absence of any bend at the warmer edge of the plateau makes uninhibited 
diffusion unlikely in depleting lithium. Moreover, an age spread of only 3 Gyrs would produce a dispersion of 0.2 dex on the 
plateau, which, again, has not been observed. 

Vauclair (1988) suggested rotation induced turbulence as the agent leading to a nuclear destruction of Li from an 
abundance as high as A(Li)=3.0 down to the observed plateau value, still approximately preserving the plateau shape. 
Similar computations with a different formulation of the rotation induced turbulence have been performed by Pinsonneault 
et al. (1992) and Chaboyer & Demarque (1994). These models predict a bending of the Li abundance on the warmer stars very 
similar to that predicted by diffusive models. In these models a considerable dispersion is expected depending on the stellar 
initial angular momenta. Only assuming an identical initial rotational velocity it is possible to deplete Li in a uniform way 
down to the observed value. Moreover, the initial rotational velocity has to be rather high, which then requires an efficient 
breaking mechanism to slow down the stars to the low rotational velocities observed. Vauclair & Charbonnel (1995) started 
from an initial value of 100 kms -1 going down to the 2 kms -1 presently observed. From an initial value of A(Li)=2.7 their 
model predicts at T c g — 6500 K a A(Li)= 1.6, with a negative slope on the plateau resembling the one obtained by pure 
diffusive models. These authors found it difficult to account for lithium dispersion less than about 20 to 50 % in the plateau 
with such models. It is also to note that an age spread of 3 Gyrs alone would produce a dispersion of 0.3 dex on the plateau. 

Li depletion induced by winds has been recently proposed by Vauclair & Charbonnel (1995). A mass loss flux between 
10" 13 and 10~ 12 MOyr^ 1 can deplete Li producing a positive slope on the plateau as a function of the stellar mass (i.e. T e g ). 
Higher mass loss rates would destroy Li entirely, while smaller rates would deplete less Li. A stellar wind of 10~ 13 M0yr~ x 
produces an almost unnoticeable Li depletion, while a wind of 10 -12 MOyr -1 produces a slope of 0.4 dex along « 1000 K; 
that is a slope of the same order of that detected by Thorburn (1994) or Ryan et al. (1996). As a comparison the solar wind 
is of 10~ 14 MQyr -1 , but there are no measurements of winds in Pop II dwarfs. The absence of such a slope as well as the 
absence of any detectable dispersion is arguing against the mass loss model. If the winds were of different rates in different 
stars they would produce a dispersion on the Li abundances, and an age spread of 3 Gyr would produce a dispersion of 0.25 
dex. 

In summary, all depletion mechanisms predict features which are not observed in the present data. The absence of a 
downturn in Li abundance at the hottest edge of the plateau and the absence of dispersion on the plateau itself argue strongly 
against significant depletion by diffusion or rotational mixing (Vauclair 1988, Pinsonneault et al. 1992); the absence of a 
significant slope with T c g and the absence of intrinsic dispersion rule out stellar winds as a possible source for Li depletion 
(Vauclair & Charbonnel 1995). Thus we consider that our results strongly support the view that the observed A(Li) in Pop 
II stars coincides with the primordial value. 

5.2 Primordial Li 

The absence of dispersion and of significant trends in our analysis, in particular when the small trend with T c g , expected 
from standard models, is considered, allows to take the mean value of Li of the 41 stars on the plateau as the value of the 
primordial Li. The weighted mean, in which each abundance is weighted inversely by its own variance, is A(Li) = 2.20 with an 
uncertainty of the mean value of only ±0.012. When we take the values corrected for the mild depletion predicted by standard 
models and NLTE effects the weighted mean is: 

A(Li) p = 2.238± 0.012 
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Figure 9. The theoretical SBBN 7 Li yields as a function of r/, following the parametrization of Sarkar (1996), are compared to our 
estimate of primordial Li. The thick solid line corresponds to A(Li)=2.238, the thin solid lines correspond to 1<t errors, while the dotted 
lines correspond to the Icr plus systematic error. The thick solid curve is the prediction of SBBN, thin solid curves are 1<t errors related 
to uncertainties in the nuclear reaction rates. 

This value is in excellent agreement with the NLTE corrected value of A(Li)=2.224±0.013 by Molaro et al. (1995a). It is 
significant that the two samples, which make use of a spectroscopically determined T e ff and of the semi-direct method of the 
IRFM, give the same average Li abundance. 

Any offset in the zero point of the different T e ff scales would affect the absolute value of the plateau, while different 
slopes in the T c g scales may introduce (or suppress) trends with T e g . The Alonso et al. (1996a) and Fuhrmann et al. (1994) 
scales are on average hotter than photometric scales and this explains why old mean values for the Li on the plateau were 
close to A(Li)=2.08 (Molaro 1991). As already mentioned in section 2, the zero-point of T e £f scales in the low main sequence 
is difficult to establish. The Alonso et al. (1996a) scale is able to reproduce both the Procyon and Sun T c s , within the 
measurement errors, and therefore no significant offsets are expected. Alonso et al. (1996a) estimate in their T ff a systematic 
error of « 1.25% at T e g = 6000 K, related to the uncertainty of the calibrations of J,H and K magnitudes. This implies a 
systematic error for Li of ± 0.05 dex, which we have considered separately since it is not eliminated in the averaging process. 
The determination of the value of primordial Li is therefore dominated by systematic rather than random errors. Still larger 
systematic errors can be hidden in our ultimate ability of modelling stellar atmospheres (Kurucz 1995). Our estimate for the 
primordial lithium is: 

A(Li) p = 2.238±0.012i CT ±0.05 sys 

In Fig. 9 the SBBN theoretical Li yields are shown, computed by using the parametrization provided by Sarkar (1996). 
The Li Pop II value corresponds to two possible values for r\ = ns/n^: rjio = l-7l '3 or r lio = 4.0iQg, where r/io = 10 10 r;. 
The minimum of the theoretical Li yield is excluded when one considers la errors in the theoretical predictions and la errors, 
plus the full systematic error, in the Li abundances. However, considering the 2a errors, both in the theoretical yields and 
in the observations (2a s tat + la s yst), the minimum of the Li curve is allowed, increasing considerably the allowed rj range to 
1.3 < nio < 5.4 . 

Our low r\ intercept coincides with the r/io = 1.8 ± 0.3 obtained from high D/H (?s 10~ 4 ) values (Songaila et al. 1994, 
Carswell et al. 1996, Wampler et al. 1996, Rugers & Hogan 1996) and is in good agreement with the rj derived from the 
primordial He, Y p = 0.228 ± 0.005 (Pagel et al. 1992), with three neutrino flavours. The high rj intercept is more consistent 
with the rjio ~ 3 implied by the primordial helium Yp=0.241 ±0.003 obtained by Izotov, Thuan & Lipovetsky (1994) by using 
revised recombination coefficients for the He I lines. 

Our low v value is inconsistent with the values nio = 6.4^7 from the low D/H=2xl0~ 5 observed in high redshift 
absorption systems (Tytler et al. 1996, Buries & Tytler 1996; but see Wampler 1996 ). To achieve consistency between the 
low D/H value and Li a 0.5 dex of depletion in Li is required, which is not supported by the present analysis. However, 
considering the 2a errors, the total allowed window for r\ overlaps with the n window implied by the low D/H observations. 

Turner et al. (1996) use the presolar value of D+ 3 He from the recent measurement of 3 He in the local ISM of Gloeckler 
& Geiss (1996) to deduce the value of T) for different assumptions on the evolution of 3 He in low-mass stars and of metal 
ejection by massive stars. Their allowed window for r\ (2 < r/io < 7) overlaps with ours, although our low rj value is close to 
their lower limit. Comparison of figure 2 of Turner et al. (1996) with our rj range seems to suggest that 3 He is preserved but 
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not produced by low-mass stars. 

The two intercepts on the Li theoretical curve correspond to two values for the baryonic density. From the relation 
n B h 2 = 0.00366j7io we have: 

n B h 2 = 0.00621°;°°^ or Q, B h 2 = 0.01461°;°°^ 

where the Hubble parameter is Ho = 100ft kms _1 Mpc _1 . Considering that the luminous matter has a contribution to O of 
Sl LUM = 0.004 + 0.0007/*T 3/2 (Persic & Salucci 1996), we have that over the entire range of allowed Ho the baryonic density 
from SBBN remains always greater than Qlum suggesting the presence of dark baryons. 
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